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Abstract: We consider twist-L operators in the planar jV = 6 superconformal Chern- 
Simons ABJM theory. Their anomalous dimension 7^^(iV) is a function of the twist L, the 
spin N, and the dressed coupling of ABJM. We show that at next-to-leading order in the 
large spin expansion, this anomalous dimension is related to that of 7\/^ = 4 SYM twist 
operators by a simple scaling law. 
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1. Introduction 

The maximally supersymmetric A/" = 4 SYM theory is an ideal theoretical laboratory for 
the investigation of planar gauge theory. At weak coupling, all-order integrability leads 
to long range Bethe Ansatz equations which are an efficient computational tool and re- 
place the usual diagrammatical expansion [||]. At strong coupling, AdS/CFT correspon- 
dence provides a non trival reformulation in terms of type IIB superstring propagating 
on AdSs X 5^ [^. This leads to quantitative predictions which are hardly achievable by 
different means. 

In this perspective, the so-called twist operators are an important probe of various 
universal quantities characterizing the gauge dynamics. At weak coupling, they are gauge 
invariant single trace composite operators belonging to a perturbatively closed 5 1(2) sec- 
tor. They are built as a matrix product of L fundamental .AA = 4 SYM scalars ip with a cer- 
tain number N of covariant derivatives spread along the chain, Ol.jv ~ Tr(^^~-^ (p). 
The parameters L and N are called the twist and spin of the operator for clear reasons. 
The dual state at strong coupling is well identified in twist-2 and is described classically 
by a rotating folded string configuration |]. 

The analysis of twist operators is based on the study of their large spin limit, iV — > oo. 
In this regime, we can study their anomalous dimension 7L,j\r(5) as a function of the twist 
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L and, of course, the planar coupling g. The next-to-leading (NLO) expansion of 7L,j\r(5) 
has the general simple form ^ ^ ^. 



lL.Ni9) = fig) (log N + j^-{L-2) log 2) + BL{g) + ■■■ , (1.1) 

where we have neglected terms vanishing at large N ^ . The coupling dependent coeffi- 
cient of the logarithm is equal to twice the conventional QCD cusp anomaly [0] and is 
commonly referred as the scaling function. It is an universal quantity which is known 
at four loops by diagrammatical methods [ p^ , |l5p. An integral equation for f{g) has been 
proposed in [ ]16[ ] (BES). It can generate the weak coupling expansion of f{g) at any desired 
order. The strong coupling expansion of f{g) is also known [ |l7| ] (see also [|T^]). 

From a physical point of view, the cusp anomaly defines what is called the physical 
coupling which measures the (universal) intensity of soft gluon radiation [|^]. On the 
other hand, the physical meaning of the subleading constant term ^^(g) is more subtle. 
It has been partially clarified in the recent analysis of Dixon-Magnea-Sterman [ |T^ (DMS). 
From the work of DMS, we know that Bl is related to a piece of the infrared divergence 
of on-shell scattering amplitudes in 7\A = 4 SYM. The singular part of the amplitudes can 
be written schematically in the following factorized form 

Air ~ exp I — ^2 I , (1.2) 

where e is the IR dimensional regulator and is the n-th times iterated logarithmic 

integral. The new function G{g) is the coUinear anomalous dimension [21, 22]. The impor- 



tant remark of DMS is that (at least for twist-2) G is the sum of a universal contribution 
Geik and the subleading constant ^2(5). The universal part is an eikonal contribution 
which can be extracted from the soft logarithms of the Drell-Yan process or from a suit- 
able rectangular light-like Wilson loop. The non-universal piece ^2(5) comes from the 
virtual contribution to the Altarelli-Parisi splitting kernel [23]. In Mellin space, this is 
precisely the subleading constant term appearing in Eq. ( ff.ip . We shall call it the virtual 
scaling function. The DMS proposal has been recently confirmed at strong coupling in [^]. 



For the case of .A/" = 4 SYM, an integral equation for Bi,{g) has been analyzed in [25, 



46]. The equation is derived neglecting wrapping effects. These are well known at leading 
order for s[(2) twist operators [|26[]. At weak coupling, they do not affect the constant 
Si (5). Also at strong coupling, this turns out to be true as an outcome of the analysis 
of[^]. 

Recently, twist operators have also been introduced and studied in the ABJM the- 
ory [ ^ which has an integrable structure quite close to that oiU = ^ SYM, despite being 
a very different theory. ABJM is a 3d U{N) x U{N) Chern-Simons gauge theory with 
opposite levels +k, —k. It has M = 6 superconformal symmetry and describes the low en- 
ergy limit of N parallel M2-branes at a C^/Zfc singular point. For large N, k and fixed ratio 



and partially investigated in string theory[ 1^ 



^The subleading corrections to Eq. (1.1) are very interesting and are related to a generalized Gribov- 
Lipatov reciprocity as first discussed in 1^, lo[ 0. It is a property very well tested at weak coupling [p^. 
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A = N/k, we identify Zk — and M theory on the gravity dual orbifold AdS^ x 
reduces to type IIA string on AdS^ x CP^. The classical string integrability ^ ^ has 
a gauge theory quantum counterpart first analyzed in [ pT| , |3^ ]. In [p3|], a set of all-loop 
Bethe-Ansatz equations have been proposed for the full osp(2, 2|6) theory consistent with 
the string algebraic curve at strong coupling [p^. The equations depend on a dressed 
coupling h{X) which takes into account the fact that the one-magnon dispersion relation 
is not protected by supersymmetry ^ ^ ^]. The conjectured 5-matrix is worked 



out m 

ABJM twist operators and their anomalous dimensions are discussed in [||. IqUI]. 



The similitudes with M = ^ SYM are many. At strong coupling, the dual string state is a 
folded string rotating in AdS^^ with large spin N and with angular momentum J ~ log N 
in CP^ [42] not so different than the AdS^ x case. At weak coupling, their anomalous 



dimension is captured by all-loop Bethe Ansatz equation which are precisely those of 
jV = 4 SYM apart from a phase twist. This phase deformation [|^] changes drastically 
some known fine features which are shared by JV = 4 SYM and QCD (Low-Burnett- 
Kroll wisdom, Gribov-Lipatov reciprocity). Nevertheless, the scaling function of ABJM is 
simply halved, apart from the unavoidable replacement g — > h{g), and we have [^] (see 
also 

f0^^h) = \f^=\h), (1.3) 

where, here and in the following, CS stands for Chern-Simons and identifies the ABJM 
case. This factor 1/2 can be understood in terms of the mode number change due to the 
phase twist in the (thermodynamical) large spin limit. 

One immediately asks the following questions. What about the subleading constants 
appearing in Eq. ( |1.1| ) ? How do they change in ABJM ? Is there a simple relation with the 
A/" = 4 SYM values ? The result of this paper is that at next-to-leading order, i.e. at the 
level of the expansion Eq. ( |1.1[ ), we can write 

= \ N), (NLO). (1.4) 

In other words, we predict 

l2%ih) = I f^=Hh) (log (2 iV) + 7B - 2 (L - 1) log 2) + ^ B^^-\h) + ■■■. (1.5) 

This conclusion is derived under the same hypothesis about wrapping we assumed in 
^^ = A SYM 2. Notice that Eq. (Q is false beyond NLO order. 

The plan of the paper is as follows. In Section^ we briefly recall a few necessary basic 
definitions. In Section ^ we present a new closed formula for the (asymptotic) anomalous 
dimension of ABJM twist-2 operators at fourth order, i.e. eight loops. In Section ^ we 
compute and collect our results for the NLO large spin expansion of the ABJM anomalous 
dimensions. In Section ^ we briefly recall the form of the NLO BES equation in JV = 4 



first positive test of this assumption has been presented in | |4l| ] for twist-1 and twist-2 operators at weak 
coupHng and leading wrapping order. 
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SYM as a necessary preliminary steps for the introduction of its modified version valid 
in ABJM. Section ^ illustrates the analysis of the one-loop ABJM Baxter equation which 
describe twist operators. The analysis is aimed at providing arguments for the scaling 
behaviour of the hole solutions to the Bethe Ansatz equations. In Section^ we report the 
jV = 4 and ABJM NLO BES equations in unified form and prove the scaling relation as a 
simple consequence. Finally, Sectionals devoted to some final comments. 

2. Twist operators in ABJM and their exact anomalous dimensions 



The all-loop Bethe equations for ABJM has been proposed in [p3[]. They are associated 
with the osp(2, 2|6) Dynkin diagram (in the fermionic rj = —1 grading) 




U4 



(2.1) 



Twist operators in the s[(2) sector are obtained by exciting the same number N of Ui and 
roots. As in the J\f = A case, we shall refer to the integer L as the twist of the operator. 
More details can be found in [4C]. 

Bethe Ansatz equations are written in terms of the deformed spectral parameters 
x'^lu) defined by 

_ 1 
~ h 

where h{X) is the interpolating coupling entering the one-magnon dispersion relation. For 
twist L operators they are 



u± 



(2.2) 



N 

n 



Uk 



Uj + I 



Uk 



„2 



(2.3) 



The only difference compared with J\f = ^ SYM is the extra minus sign whose effects are 
discussed in [|^, The factor ctbes is the Beisert-Eden-Staudacher dressing phase. The 
momentum constraint is automatically satisfied for Bethe root distributions symmetric 
under u —uas those we are interested in. 

The contribution to the energy /anomalous dimension from the Asymptotic Bethe 
Ansatz (ABA) equations is conveniently written in terms of 



p(u) = —i log 



x^{u) 
x~{u) ' 



^ cot^ \/l + 16 sin^l, 



(2.4) 



and reads 



N 



k=l 



1 + 16 /i2 sin^ — - 1 
2 



(2.5) 
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From the results of [41], we know the closed form of the asymptotic anomalous dimension 
at three orders for both twist 1 and 2. We recall here the expressions for the reader's 
benefit. With the usual definition of (nested) harmonic sums 



«.w = E^^. ^....m = f:Ms,.(„). (2.6) 



we have for twist 1 



Ti^fW =4(5i-5_i), (2.7) 

- 5_i,_i,i - 5i,_i,_i + 5i,_i,i), (2.8) 



and a much longer expression for JieiN) which can be found in [41]. All harmonic sums 
have in this case argument 5a = Sa{N). For twist-2, we find instead the compact expres- 
sions 



72°i(^) =4 5i, (2.9) 

j^^^iN) = 4 53 - 8 5i,2 - 4 52,1, (2.10) 
l2fiW = 8 55 - 24 5i,4 - 32 52,3 - 20 53,2 - 16 54,1 + 32 5i,i,3 + 24 5i,2,2 + 

+24 5i,3,i + 20 52,1,2 + 24 52,2,1 + 8 53,1,1 - 16 5i,i,2,i. (2.11) 



and this time the argument is = Sa{N/2). 

These expressions do not include wrapping corrections. We shall compute the virtual 
scaling function from these asymptotic expressions assuming, as in .A/" = 4 SYM that it is 
independent on wrapping. 



3. The fourth order twist-2 anomalous dimension 



As remarked in [^] and in close analogy with twist-3 fields in 7\A = 4 SYM, the twist-2 
anomalous dimensions involve nested harmonic sums with positive indices only. Such 
a restricted maximal transcendentality Ansatz reduces a lot the complexity of the deter- 
mination of higher orders. In particular, we have determined the 8-th loop asymptotic 
anomalous dimension where dressing first appear. This will be a useful ingredient for the 
check of our results for the virtual scaling function. After a straightforward computation, 
we find 

72^8'(^) = G^U^) + C3 (3.1) 
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where 



G^j = 4(5 57 - 20 5i,6 - 39 53,5 - 45 53,4 - 35 54,3 - 27 55,2 - 15 Se.i + 48 5i,i,5 + 
+62 5i,2,4 + 58 5i,3,3 + 52 5i,4,2 + 38 5i,5,i + 58 52,i,4 + 64 52,2,3 + 
+68 52,3,2 + 55 52,4,1 + 44 53,1,3 + 54 53,2,2 + 52 53,3,1 + 31 54,1,2 + 33 54,2,1 + 
+26 55,1,1 - 48 5i, 1,1,4 - 52 5i, 1,2,3 - 60 5i, 1,3,2 - 52 5i, 1,4,1 - 40 5i,2,i,3 + 
—58 5i,2,2,2 — 58 5i,2,3,i — 32 5i,3,i,2 — 42 5i,3,2,i — 32 5i,4,i,i — 32 52, 1,1,3 + 
—54 52,1,2,2 — 54 52,1,3,1 — 38 52,2, 1,2 — 52 52,2,2,1 — 42 52,3, 1,1 — 18 S'3, 1,1,2 + 
—32 53,1,2,1 — 34 53,2,1,1 — 12 54,1,1,1 + 32 5i, 1,1,2,2 + 32 5i, 1,1,3,1 + 16 5i, 1,2, 1,2 + 
+32 5i, 1,2,2,1 + 24 5i, 1,3, 1,1 + 4 5i,2, 1,1,2 + 24 5i,2, 1,2,1 + 28 5i,2,2,i,i + 
+20 52,1,1,2,1 + 28 52,1,2,1,1 + 12 52,2,1,1,1 — 16 5i, 1,1,2, 1,1), (3.2) 

^CS,dressing ^ 33 ^ _ _ 3 3, , ,). (3.3) 



4. The virtual scaling function in ABJM from weak coupling 

Up to vanishing terms for N oo, the large spin expansion of our exact anomalous 
dimensions can be computed by the tricks reported in Appendix (A). Partial results can 



be found in [l4l[]. We add here the 6 loop twist-1 and 8 loop twist-2 results. They read 



j2ii{N, h) = f'^^ih) (log + 7E + log 2) + B2l,{h) + ■■■, (4.1) 
j2i2{^, h) = f^^Hh) (log iV + 7B - log 2) + B2i2{h) + ■■■, (4.2) 

where dots denote terms vanishing for large spin. The CS scaling function is half the 
jV = 4 value 

f^^{h) = lf^=\h), (4.3) 
f^-\g) = 8 5^ - ^ TT^ 5^ + ^ TT^ / - 16 TT^ + 4 Cl) / + 

/ 887 a ,2 /2 I /in / / ^ „10 



+ 32 (^_-^7r« + -7r^CI + 40C3C5j g'' + ■ ■ ■ ■ (4.4) 
The virtual scaling function derived from harmonic sums turns out to be 



dCS 
■°L=1 



+ ( ^ + 80C5 ]h' + (h') , (4.5) 



B2i2 = 4C3/i' - 88C5/i' + (-^^'Cs + 167r2C5 + lUOCr) + O [h'") (4.6) 
These values must be compared with the known weak-coupling expansion valid for twist- 
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L fields in j\A = 4 [H, E|] 



B^=\9) = 8 (2 L - 7) Cs 5' + (- 3 - 4) TT^ Cs - 8 (21 L - 62) Cs j / + 

+ ("^ (3 L - 13) TT^s + ^ (111' - 32) 7r2C5 + 40 (46 L- 127) Cr) / + 



+ -64(2L-7)C| 



128 8 

— (11 L - 49) Cs - ^ (103 L - 310) tt^ C5+ 



40 



(25 L - 64) TT-^ C? - 392 (55 L - 146) Cq 5 " + 



(4.7) 



5. The NLO BES equation in M = A 

Let us recall how B-^~^ is computed. In the by now standard notation of the BES pa- 
per []l6|], one has to solve the integral equation for the quantity a{t) closely related to the 
Fourier transform of the Bethe root density and obeying 



ait) 



t 



e* - 1 

\ /-oo 



dt' 



K{2 g t, 0) (log + 7e - (L - 2) log 2) 
2 L-2 



e*' - 1 e*'/2 + 1 



8gH 

{K{2gt,2gt')- K{2gt,0)) + 



iM2gt)-l) + 



-452 / dt' K{2gt,2gt')a{t') 
Jo 



The anomalous dimension is simply given by 

I^=\g) = l6g'ai0). 
The kernel appearing in the above equation is 

K{t, t') = Koit, t') + Ki{t, t') + Kdit, t'), 



with 



= ^ E(2^-l)-^2n-l(^)J2n-l(<'). 



n>l 



^1 = ^ T.i^ri)J2n{t)J2n{t'), 



roo 



Kd{t,t') = 852 / dt" K,{t,2gt")-j^Koi2gt",t'). 
Jo - 1 



(5.1) 

(5.2) 

(5.3) 

(5.4) 
(5.5) 
(5.6) 



Now, the crucial point is that Eq. ( p31 ) is in close relation with the general properties of the 
transfer matrix which appear in the Baxter equation of jV = 4 as discussed in p5|]. Since 
we want to derive a suitable modification of Eq. ( |5.1[ ) valid in ABJM, we now analyze the 
Baxter equation for the ABJM twist operators. 
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6. Analysis of the ABJM sl(2) Baxter equation 



The following discussion will be done at the one-loop level will turn to be enough for our 
purposes. It is the adaptation of the methods described in p7[ ] to the ABJM (twisted) case. 
The starting point is the Baxter (5(«)-operator where u is the spectral parameter. It acts on 



the L-sites long s[(2) spin chain and obeys the (twisted) Baxter equation pSQ 

2 



u + ^] Q{u + i)-[u-^] Q{u-i) = iXu)Q{u). (6.1) 



The operator t{u), also known as the auxiliary transfer matrix, is a polynomial of degree 
L — 1 in u with coefficients being local charges g„ also acting on the spin chain states. The 
more precise form of t{u) is 

L 

F(u) = z(2iV + L)u^-^ + ^ (6.2) 

n=2 

Both the Baxter operator and the transfer matrix can be simultaneously diagonalized with 
the Hamiltonian. Replacing the operators with their eigenvalues, we obtain the same dif- 
ference equation for the Baxter function Q{u). The suitable boundary condition is simply 
the requirement that Q{u) is a polynomial of degree N 

N 

Qiu)=J^ Yl{u-Un). (6.3) 

n=l 

Its roots are readily identified with the Bethe roots since they obey the Bethe Ansatz equa- 
tions. Notice that the operator Q is required in order to compute the precise form of the 
energy eigenstates. Instead, for the eigenvalues, the polynomial Q{u) is enough. 

In practice, given a certain spin N, we can replace in the Baxter equation a polyno- 
mial Ansatz with undetermined charged g„. Solving for the polynomial coefficients, we 
end with a set of algebraic equations for the charges. Each solution is associated with 
a specific sl(2) module appearing in the subset of [—1/2]®^ states compatible with the 
twisted boundary conditions. The solutions 5„ of t{dn) = are called holes. They are dual 
solutions of the Bethe equations whose role will be clarified in a moment. 

If we look for an even Q[—u) = Q{u), then the transfer matrix (eigenvalue) has also a 
definite parity t{—u) = (— 1)^"-^ t{u). Thus, the cases L = 1, 2 are particularly simple. For 
L = 1 we do not have holes nor charges. For L = 2 we have a single null hole 5 = and, 
again, no charges. For L > 2, non-trivial charges start to appear. As discussed in [ ]47| , 
the behaviour of the spin chain in the regime iV 3> ij 3> 1 is largely determined by the 
N dependence of the holes. In particular, it is important to determine whether the holes 
grow large as i\r oo or vanish. In 7\A = 4, it is known that two holes have a size which is 
(D{N) while the other L — 2 holes are small in the sense that they vanish like 1/ log N. 

We claim that in the ABJM case, all holes are 0(1/ log N). This statement can be effi- 
ciently checked by the methods of [ p7[ ] which provide a good approximation to the holes 
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(accurate enough for our purposes) in the N ^ L ^ 1 regime. The main idea is that in 
the full region u = 0{{2N + L)°) one can solve separately the two half-Baxter equations 



U + -] Q+{u + i) = +t{u)Q+{u), (6.4) 



u--] Q-{u-i) = -t[u)Q-[u). (6.5) 



where we write 



The solution is 



L-l 



t{u) =i{2N + L) ]\{u- 5n). (6.6) 



L-l 



Q+(u) = [(2i\r + L)z]-" -— — Y{V{i5n-iu) (6.7) 

with Q- being related to <5+ by conjugation, regarding u as real. The approximate asymp- 
totic Baxter function is then simply 

Q^{u) = Q+{u) g-(-z/2) + Q_{u) Q+(z/2). (6.8) 

The request of zero residue at the holes leads to the quantization condition 

5n log(2 N + L) + LAr^v(^--i5n) + Arg r(l + i J„ - i J, ) = | (6.9) 
where kn are integers. Asymptotically this means 

^•■ = f log(^ + L/2) + f2L + l)log2 + ,. ="W'°«^'- 



We can test the accuracy of Eq. ( |6.9|) by comparing its solutions with the exact holes 
j-jexactj. jj-jg results are shown in the Tables reported in Tab. (1-2) for the non-trivial 
cases L = 3, 4, 5. The agreement is very good. We thus have very strong indications that 
all the L — 1 holes are vanishing for large spin. As an immediate (one-loop) consequence, 
we can evaluate the anomalous dimension which (from Eq. ( |6.8[ )) is 



7£|(iV) = 4 1og(2iV + L)-4./;(l) + 2 ^ 



^[l^+l5n]+^[--i5n] -2V(1) 



, (6.11) 



where the second —4 ■^(1) term comes from the different numbers of factors in the last two 
terms of Eq. ( ^^ and is absent in A/" = 4. Taking the large N limit with 5„ — ^ we find 



^2%{N) = 4 [log iV + 7E + (3 - 2 L) log 2] + ■ ■ ■ (6.12) 
In particular, this is in agreement with the explicit one-loop results at L = 1, 2. 
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7. The integral equation determining BL{g) 



From the results of the previous section, we can repeat the steps leading to Eq. The 
main differences are listed below. 

1. In TV^ = 4 SYM, the natural universal constants accompanying log N were the com- 
bination 

log i\r + 7b - (L - 2) log 2. (7.1) 
This is replaced in ABJM by the (one-loop) combination 

log iV + 7e - (2 L - 3) log 2. (7.2) 



2. The explicit factor L — 2 in Eq. ( [5.1D can be traced back to the number of small holes 
Nh which vanish for large N. It becomes L — 1 in ABJM. Similarly, there are various 
terms with factors 2 which are really L — Nh- This is 1 in ABJM. 



In conclusion, we can write a modified Eq. ( p.l| ) depending on a parameter ^ such that 



e^=^ = l, ^' = 1- (7.3) 



It reads 



2 



iK{2gt,Q) (^log|^ + 7B- (^|-2) log 2^ - ^ ( Jo(2 5 1) - 1) + 
L [-^l-^m^l] i^i^9t,2gt')-K{2gt,0)) + 



00 



/ dt' K{2gt,2gt')a^{t') 



(7.4) 



From this equation we immediately obtain the following exact scaling relation valid at 
NLO in the large spin expansion 

1lHn)=^-^^\2N). (7.5) 
As a corollary, we obtain Eq. (^) which is indeed verified by our explicit data at L = 1, 2. 

8. Comments 

In this paper, we have proved the scaling relation Eq. (pT^. This is not a totally surprising 
result. Indeed, it already appeared in the discussion of a rigid circular string stretched in 
both AdSi and along an of CP^ and carrying two spins [^]. It has also a counterpart 
in the case of the generalized scaling function of [iS] as remarked in p3|]. Here, we have 



given a proof valid in the case of twist operators belonging to a special s[(2) subsector 
with many similarities to the M = ^ SYM case. We have studied in details the necessary 
changes in the derivation of the integral equation for . This required the analysis of 
the asymptotic properties of the hole solutions to the twisted Baxter equation describing 
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this sector of ABJM. Eq. ( pT^ is checked against the exact known anomalous dimensions 
of twist operators in ABJM presented in as well as the new fourth order result for 
twist-2 illustrated here. Notice also that Eq. ( |7.5D is valid at all orders in the coupling. In 
particular it permits to obtain the strong coupling expansion of the virtual scaling function 



from the recent results of [25] with no effort. 

As a final comment, let we remark that Eq. ( |L^ ) answers a technical problem, i.e. the 
role of the phase deformation of ABJM at the level of the virtual scaling function. This is 
certainly interesting, but deserves a more sound physical motivation. We mentioned in 
the Introduction the important recent developments connecting the virtual scaling func- 
tion of 7^ = 4 SYM to the properties of the on-shell scattering amplitudes. This connection 
has not been explored in ABJM. Compared to jV = 4 SYM, the quite different nature of the 
gauge theory side sets itself against the close integrability structure and, in our opinion, 
makes the problem an intriguing one. 
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A. Next-to-leading large spin expansion of harmonic sums 

Let us remind a simple trick which reduces the NLO expansion of harmonic sums to the 
evaluation of multiple values. The general large N expansion of a nested sum Sa{N) has 
a singular part which is a polynomial in log N, a constant term, and a remainder which 
vanishes at iV = oo. The degree of the singular polynomial equals the number of leading 1 
indices in a = (ai , . . . , a^). This singular part can be extracted by using the shuffle algebra 
described in [||]. The main point is that if we start with (ai 7^ 1, but it can also be absent) 

Si,l^^^,ai,a2,...,, (A.l) 

k 

we can add and subtract the product 

-^l 5i,...,i,ai,a2,...,- (A.2) 



k + 1 



Using the shuffle algebra in one of this terms, we cancel the initial sum in Eq. (Q. Re 



peating iteratively this procedure we obtain the desired result. As an example let us con- 
sider 5i, 1,1,2. Applying the above algorithm, we prove that 

'^1,1,1,2 = 7'52'S'l + ^) Sf + - 53,1 + 52,l,l) Si + 



6 ' V 2 2 J ' V 2 
6 3 6 2 ■ 2 



02,3 03,2 04,1 ^ a n q /A a^ 

+— H — s ;:; — r -^2,1,2 + -^3,1,1 - ^2,1,1,1. K^-^) 
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Taking the values at infinity we simply get 

Si 1 1 2{N) = — TT^L^ - -CsL'^ + — TT^L - C5 - — TT^Cs + " " " , (A.4) 

where L = log iV + 7e is the NLO expansion of Si{N). 
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N 


'5l=3 




<^L=4 




4 


0.22795472 


0.22787539 


0.40848291 


r\ A r\ c\ r\ a r> 

0.40824222 


8 


0.21022561 


0.2102162 


0.37840139 


0.37832836 


12 


0.20032517 


0.20032661 


0.3613073 


0.36127663 


16 


0.1936300 


0.19363366 


0.3496743 


0.34965917 


on 


U.iooo4//i 




u.o4uyyu/ 


U.o4Uyo/c)4 


24 


0.18471752 


0.1847212 


0.33413129 


0.33412668 


28 


0.18149529 


0.18149861 


0.32850149 


0.32849884 


32 


0.17877816 


0.17878113 


0.32375148 


0.32374999 


36 


0.17643811 


0.17644076 


0.31965923 


0.31965845 


40 


0.17438935 


0.17439172 


0.31607559 


0.31607526 



Table 1: Comparison between the exact holes and the solutions to the quantization condition (qc) 

for 1/ = 3, 4. In the first case, we have two opposite holes ±5. In the second, we have three holes 
0, ±5. We show the non trivial value 5 in all cases. The integers fc„ entering the quantization 
condition are ±1 for L=3, and ±2, for L = 4. 
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rex 


.qc 






4 


0.15014187 


0.15013878 


0.57240546 


0.57206928 


8 


0.14425421 


0.1M2A979 


0.53079427 


0.53066078 


12 


0.14046208 


0.14045861 


0.50680642 


0.50673907 


16 


0.13769667 


0.13769405 


0.49041072 


0.49037163 


20 


0.13553582 


0.13553381 


0.47815412 


0.4781293 


24 


0.13377116 


0.13376958 


0.46846905 


0.46845227 


28 


0.13228511 


0.13228384 


0.46052135 


0.46050948 


32 


0.13100508 


0.13100404 


0.45381825 


0.45380955 


36 


0.12988321 


0.12988235 


0.44804621 


0.44803966 


40 


0.12888638 


0.12888565 


0.44299429 


0.44298925 



Table 2: Comparison between the exact holes and the solutions to the quantization condition (qc) 
for L = 5. The holes are ±5l^s,±{5l^5)' . We show the positive holes. The integers fc„ entering 
the quantization condition are ±3, ±1. 
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